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Abstract 

We present a general gauge invariant formalism for defining cosmological averages that 
are relevant for observations based on light-like signals. Such averages involve either null 
hypersurfaces corresponding to a family of past light-cones or compact surfaces given by 
their intersection with timelike hypersurfaces. Generalized Buchert-Ehlers commutation 
rules for derivatives of these light-cone averages are given. After introducing some 
adapted "geodesic light-cone" coordinates, we give explicit expressions for averaging the 
rcdshift to luminosity-distance relation and the so-called "redshift drift" in a generic 
inhomogeneous Universe. 



1 Introduction 



It is by now well-known (see, for example, [1]) that averaging solutions of the full inho- 
mogeneous Einstein equations leads, in general, to different results from those obtained 
by solving the averaged (i.e. homogeneous) Einstein equations. In particular, the averag- 
ing procedure does not commute with the non-linear differential operators appearing in the 
Einstein equations and, as a result, the dynamics of the averaged geometry is affected by so- 
called "backreaction" terms, originating from the inhomogeneities present in the geometry 
and in the various components of the cosmic fluid. 

Following the discovery of cosmic acceleration on large scales, interest in the possible 
effects of inhomogeneities for interpreting the data themselves has considerably risen (see 
[2, 3] for some recent review papers). It has even been suggested (see e.g. [4]) that the 
dynamical effects of the backreaction could replace, at least in part, the dark-energy sources 
as an explanation of such a cosmic acceleration, thereby providing an elegant solution to 
the well-known "coincidence problem" . 

As a consequence, much work has been done in these last few years on trying to formulate 
a suitable "averaged" description of inhomogeneous cosmology. In most of these works, 
following Buchert's seminal papers [5], the effective geometry emerging after the smoothing- 
out of local inhomogeneities has been determined by integrating over three-dimensional 
spacelike hyper surfaces. 

However, as pointed out long ago [6, 7], a phenomenological reconstruction of the space- 
time metric and of its dynamic evolution on a cosmological scale is necessarily based on 
past light-cone observations, since most of the relevant signals travel with the speed of light. 
Hence, the averaging procedure should be possibly referred to a null hypersurface coinciding 
with a past light-cone rather than to some fixed-time spacelike hypersurface. Nonetheless, 
such a light-cone averaging procedure, whose importance has been repeatedly stressed in 
the specialized literature (see e.g. [2, 8]), has never been implemented in practice. 

The aim of this paper is to introduce a general (covariant and gauge invariant) prescrip- 
tion for averaging scalar objects on null hypersurfaces, to apply it to the past light-cone of 
a generic observer in the context of an inhomogeneous cosmological metric, and to provide 
the analog of the Buchert-Ehlers commutation rules [9] for the derivatives of light-cone av- 
eraged quantities. We will also introduce an adapted system of coordinates (defining what 
we call a "geodesic light-cone frame" , which can be seen as a particular specification of the 
"observational coordinates" introduced in [6, 10]) where the averaging prescription greatly 
simplifies, while keeping all the required degrees of freedom for applications to general 
inhomogeneous metric backgrounds. 
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In order to illustrate our light-cone averaging procedure we will propose, in particular, 
two possible physical applications: the averaging of the luminosity-redshift relation and 
that of the so-called "redshift drift" parameter [11] (see [12, 13] for recent discussions), 
both evaluated for a generic inhomogeneous cosmological geometry. We will concentrate 
our attention on the realistic case in which the reference observer is geodesic, having in 
mind a model of light ray propagation based on the geometrical optics approximation (see 
e.g. [14, 15]). 

Finally, we stress that our intent, in the present paper, is not to give a detailed physical 
discussion and interpretation of the many phenomenological effects of averaging, and to 
bring arguments supporting (or disproving) plausible theoretical alternatives to the "stan- 
dard" dark-energy /cosmic acceleration scenario. We simply present a first step towards such 
an ambitious program, providing a formal procedure allowing an automatic implementation 
of light-cone averaging, a procedure that - to the best of our knowledge - was missing in 
the present literature, and that we believe to be of some utility because of its covariance 
and model independence. A discussion of the possible phenomenological consequences of 
its application to specific cosmological scenarios will be presented in a forthcoming paper. 

The paper is organized as follows. In Sect. 2 we recall the averaging prescription intro- 
duced in [16, 17] and we generalize it to different kinds of averaging on null hypersurfaces. 
In Sect. 3 we give the covariant and gauge invariant version of the Buchert-Ehlers commu- 
tation rules for light-cone averages. In Sect. 4 we introduce a coordinate system in which 
the light-cone averaging prescriptions and commutation rules take a much simpler form. In 
Sect. 5 we discuss a possible approach to the average of the luminosity-redshift relation and 
of the redshift drift on the light-cone of a geodesic observer. Our conclusive remarks are 
presented in Sect. 6. 



2 Gauge invariant light-cone averaging 

Let us first briefly recall the approach given in [16, 17] to gauge invariant averaging on a 
three-dimensional spacelike hypersurface S(-A), embedded in our four-dimensional space- 
time Assuming the hypersurface (or a spacelike foliation) to be defined by an equation 
involving a scalar field with timelike gradient A(x): 

A(x)-A = 0, (2.1) 



2 



the gauge (and hypersurface-parametrization) invariant definition of the integrai of an ar- 
bitrary scalar S(x) and of its average on such hypersurface was given in [16, 17] as 1 : 

I(S; A ) = f d A xJ-g{x) S(A(x) - A ) J-d„Ad^A S(x); 
J Ma 

{S)Ao = M (2 - 2) 

(we are using the signature (—,+,+,+)). Here the spatial hypersurface has no bound- 
ary. However, as shown in [16, 17], a possible spatial boundary can be added through the 
following extension of the previous integral: 

I(S;A ;B )= [ d 4 xJ -g(x)5(A(x) - A )@(B - B{x))J -d^A S(x) , (2.3) 

J Mi 

and similarly for the corresponding average (6 is the Heaviside step function, and B is a 
positive function of the coordinates, with spacelike gradient). As discussed in [16, 19], this 
is still a gauge invariant expression if B(x) transforms as a scalar, while it gives violations 
of gauge invariance if B is not a scalar and keeps the same form in different coordinate 
systems. Even in that case, however, gauge invariance violations go to zero when we choose 
Bq in such a way that the size of the spatial region goes to infinity [16, 19]. 

The above procedure unfortunately fails if A(x) = Aq defines a null (lightlike) hyper- 
surface, since in that case d^Ad^A = 0. In order to circumvent this problem let us start 
with a spacetime integral where the four-dimensional integration region is bounded by two 
hypersurfaces, one spacelike and the other one null (corresponding e.g. to the past light- 
cone of some observer). Let us choose, in particular, the region inside the past light-cone of 
the observer bounded in the past by the hypersurface A(x) = Aq: clearly a gauge invariant 
definition of the integral of a scalar S(x) over such a hypervolume can be written (in a 
useful notation generalizing the one used above) as 

I(S;-,Aq,V )= I d 4 x^g-e(V -V)@(A-A )S(x), (2.4) 

J Mi 

where V(x) is a (generalized advanced-time) scalar satisfying d^Vd^V = 0, where Vq specify 
the past light cone of a given observer, and where the "— " symbol on the l.h.s. denotes the 
absence of delta- like window functions. 

Starting with this hypervolume integral we can construct gauge invariant hypersurface 
and surface integrals by applying to it appropriate differential operators - or, equivalently, 
by applying Gauss's theorem to the volume integral of a covariant divergence. An example 

1 In [17] the prescription introduced in [16] is used to give a covariant and gauge invariant generalization 
of the effective equations presented in [5]. Such a generalization has been recently used to deal with the 
backreaction of quantum fluctuations in an inflationary model [18]. 
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of the latter, if we are interested into the variations of the volume averages along the flow 
lines normal to the reference hypersurface S(-A), is obtained by replacing the scalar S with 
the divergence of the unit normal to S, 

QilA n/ = -l, (2.5) 



V-d u Ad»A' 
and leads to the identity: 

/ d^x^-g e(V - V)@(A - A))V^v = 

JMa 

= - f d 4 x^@(V - V)S(A - A^J-d^AdtA 

JMa V 

+ / d 4 xV^S(V - V)@(A - Ao) * (2.6) 

JMa yf-O u Ad u A 

Hence, if we start from Eq. (2.4), and we consider the variation of the average integral by 
shifting the light-cone V = V$ along the flow lines defined by n^, we are led to define the 
hypersurface integral (with positive measure): 

I(1;V ;A )= [ d i x^-gS(V -V)e(A-A )^2^. (2.7) 
JMa V-duAd^A 

Similarly, if we consider the variation of the average integral by shifting the hypersurface 
A = Aq (along the same flow lines defined by n M ), we are led to another hypersurface 
integral: 

I(1;A ;V )= [ d A x^g~Q(V - V)5(A - Ao)J-dpAd»A. (2.8) 

J Ma 

In the first case, Eq. (2.7), the integration region is on the light-cone itself, and it is spanned 
by the variation of £(A)) along its normal, at fixed light-cone Vq (see Fig. 1, (a)). In the 
second case of Eq. (2.8) - which gives exactly the same integral as in Eq. (2.3) with V 
replacing B - the hypersurface S(^4o) is kept fixed, and the integration region describes the 
causally connected section of S spanned by the variation of the light-cone hypersurface (see 
Fig. 1, (b)). 

Further differentiation also lead to the following invariant surface integral 

I(l;V ,A Q ;-)= [ d 4 x^ 6(V - V)S(A - A )\d t yd»A\ , (2.9) 

JMa 

with a compact, two-dimensional integration region defined by the intersection of £(A)) with 
the light-cone Vq (Fig. 1, (c)). This integral, as well as the integrals of Eqs. (2.7), (2.8), is 
not only covariant and gauge invariant but also invariant under separate reparametrizations 
of the scalar fields A — > A(A) and V —> V(V). Eq. (2.9), in addition, is a particular case 
of an invariant integration over an arbitrary codimension-2 hypersurface defined by the 
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truncated light cone causally connected sphere 




(c) I(l;V ,A ; ) 

Figure 1: Three different light-cone averaging prescriptions, (a): the average of Eq. (2.7) 
on the light-cone itself starting from a given hypersurface in the past; (b): the average of 
Eq. (2.8) on the section of the hypersurface A(x) = Aq which is causally connected with 
us; (c): the average of Eq. (2.9) on a 2-sphere embedded in the light-cone. 

conditions A^ n \x) = 0, n = 1,2. In general, and in D spacetime dimensions, such an 
integral can be written as 

[ d D x^]j5( y A^\x))J\detgP m \ ; g pm = d^A^d v A^g^ , (2.10) 

(as can be shown by considering the induced metric on the (D — 2)-hypersurface), and is 
invariant under the more general reparametrizations 

A (i) i(i) (^w^ 2 )) and A^ -)■ 
^(2) (A^\A^y It can be easily checked that Eq. (2.10) reduces to (2.9) if D = 4 and 
if A^ l > = A - A and A^ = V - V are scalar functions with timelike and null gradient, 
respectively. 

In order to make contact with Eqs. (2.7), (2.8), it may be useful to remark that the 



5 



integral (2.9) can be also obtained starting from the hypervolume integral (2.4) by consid- 
ering the variation of the volume average along the flow lines normal to T>(A) for both 0(A) 
and &(V), namely by using the following window function: 

-n»V^Q(A(x) - A)KV M 9(Vb - V{x)) = 

= ^-d»Ad»A5{A(x) - A ) j_* d MltA W - V(x)) . (2.11) 

We note, finally, that averages of a scalar S over different (hyper) surfaces are trivially 
defined - with self explanatory notation - by: 

{S)VoAo = K1;V ,A ;-) ' (2 - 12) 

{S)v ° ~ mv -A ) ' (2 - 13) 

Vo 7(5; A ;V ) 

{S)A ° ~ I(1;A ;V )- (2 - 14) 

3 Buchert-Ehlers commutation rules on the light-cone 

For the phenomenological applications of the following sections we are interested, in partic- 
ular, in the derivatives of the averages defined in (2.12) and (2.13). To this purpose, let us 
first consider the derivatives of I(S; Vo,Ao; — ) and I(S; Vq;Aq) with respect to Ao and Vo 
(quantities that, like the starting integrals themselves, are covariant and gauge invariant). 

We will use the identities (here k^ = d^V): 

S(V _ V ^^AB^M , = (3,) 

and the relation k^d^Sty — Vo) = k^k^d^V — Vo) = 0. Integrating by parts, we then obtain 
d w „ „ . , fk-dS „ . \ / V-k 



aAo I(S; V , A ; -) - / ( ; V„, A a; -) + I (^±S; V , A„; -) 
^ -I(S; Vo; Ao) = I V ; A j + I 



BAo 



vk _lk»d,({dA?) 



2 (dAf 



k^A^ A »)> 
(3.3) 



where, to simplify notations, we have introduced the following definitions: k^d^S = k ■ dS, 
k^d^A = k-dA, d^Ad^S = dA ■ dS, d^Ad^A = (OA) 2 , = V • k and □ = WV^. 
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Using these results we can write a generalized version of the Buchert-Ehlers commutation 
rule [9] for light-cone average. In the case of averages defined over the integration domain 
(2.9) we find 



d 
d~A ( 



( s )v ,a 



k-dS \ 
k-dA/ 



+ 



V ,Ao 



V-k 
k-dA' 



VoAo 



V-k\ 
k-8A/ 



VoAo 



(S) VoAo , (3.4) 



while for 3-dimensional averages over the domain (2.7) we obtain: 



dA 



k-dS \ A ° 

k-dA/ 



Vb 



V-fc- 



1 k»d» ({dAf) 
V ' k 2 {dAf 



k-dA, 



A 



1 k^df, {{dAf 

2 {dAf 



k-dA 



(S) 



A 
Vb • 



(3.5) 



Vb 



Following a similar procedure we can also evaluate the generalization of the Buchert- 
Ehlers commutation rule for the derivative of I{S; Vb, Aq\ — ) and I{S; Vb; Aq) with respect 
to Vq. As we shall see in Sect. 5, such derivatives, together with the previous ones, may- 
be useful in evaluating an averaged version of the "redshift drift" parameter. After some 
calculations we are led to 



d 



( S )v ,A 



+ 



dA-dS \ 

k-dA / VojAo 

OA - V u [ fc" 



aA-Vj fe" 



[k-dS 

{dAf ) 
k-dA) 

(dAf ) 
k-dA 



{dAf 

{k ■ dAf 

s " 
k-dA l 

1 y 



V ,Ao 



V ,Ao 



k-dA, 



(S) 



V ,A > 



(3.6) 



V ,Ao 



and 



d 
dVo 



9A-dS \ A ° _ I {dAf 

k-dA ? Vo y 0b {k-dAf 



aA- d^Ad^ln^dAf^j 



V-k 



{k-dAf + 2 K ^ 



Vb 

S \ A ° 
k-dA/ Vo 

(dAf ) 

{k ■ dAf ) 



A 



s 



Vb 



+ 



OA — d^Adu In ((dA) 
, (dAf 



1 \ A ° 

k^aXU {S)v ° 



{dAf 
{k ■ dAf 



A 
Vb 

A 



(S) 



Vb 



(3.7) 



Vb 



In a similar way one could also derive equations for the derivatives of the average (S 1 )^ . 
However, since we will not discuss applications of those, we will not give their explicit form 
here. 
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4 Geodesic light-cone coordinates 



4.1 Definition of geodesic light-cone gauge 



We now turn to a special (adapted) coordinate system, = (w,T,6 a ), a = 1,2, in which 
the previous equations take a simpler form. In this sense they will play a role similar to the 
one played by synchronous gauge coordinates for spatial averages with respect to geodesic 
observers [19]. We are interested in coordinates such that the level sets of one of them define 
the past light-cones, while those of another coordinate define a set of geodesic observers. 
We claim that such coordinates, that we will call geodesic light-cone (GLC) coordinates, 
are defined by the metric: 



ds 2 = T 2 dw 2 - 2T dwdr + -f ah {de a - U a dw){d9 b - U b dw) ; 



a, b =1,2. (4.1) 



This metric depends on six arbitrary functions (T, the two-dimensional vector U a and the 
symmetric tensor 7^) and corresponds to a complete gauge fixing (modulo residual transfor- 
mations involving non-generic functions of all the coordinates) of the so-called observational 
coordinates 2 discussed in detail in [6, 10]. In matrix form, the metric and its inverse are 
given by: 



9»v 



( T 2 + U 2 
-T 

V -U a 



-r -u b \ 


7ab / 



-1/T 
-1/T -1 
-U a /T 



\ 

-U h /T 



7 



ab 



(4.2) 



/ 



where 7 ab is the inverse of 7 a {,. Clearly w is a null coordinate (i.e. d^wd^w = 0). More 
interestingly, one can check that d^r defines a geodesic flow, i.e. that 



(^r)V, (0„t) = O, 



(4.3) 



as a consequence of g TT = — 1. Thus an observer defined by constant r spacelike hypersur- 
faces is in geodesic motion. Also note that ^J z^ g = T^/pyf, with 7 = det7 a b. 

In order to understand the geometric meaning of these variables, it is useful to consider 
the limiting case of a spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) Uni- 
verse, in the cosmic time gauge, with scale factor a(t). Such a limit is easily reproduced by 
Eq. (4.1) by setting 



w = r + rj 1 t = t, T = a(t), 

Jabd6 a d6 b = a 2 (t)r 2 (d9 2 + sin 2 9d(p 2 ), 



U a = 0, 



(4.4) 



2 Note that our coordinates 9 a are not necessarily "observational", in general, but they can be reduced 
to such form (e.g. to standard spherical coordinates, parallelly propagated along the observer world-line) 
through an appropriate relabelling of null generators [10] . 
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where r] is the conformal time of the homogeneous metric: drj = dt/a. 

In these coordinates a past light-cone hypersurface is specified by the condition w = 
wq, and the momentum of a photon traveling on it, being proportional to k^ = d^w, is 
orthogonal both to itself (k^k^ = 0) and to the 2-surface generated by d$a [d^k^ = 0). 
The velocity of a generic observer defined by a scalar A as in Eq. (2.5) satisfies: 

^"=k^m- (4 - 5) 

For a geodesic observer A depends only on r and we can always set A = r, so the above 
relation simplifies to: 

k^ = T" 1 . (4.6) 

These equations will be used in Sect. 5 in connection with the redshift and luminosity 
distance. We also note for later use that, in these coordinates, the covariant divergence of 
k^ takes the simple form: 

V^ = ~d T (hlyfw) . (4.7) 

4.2 Average equations in the GLC gauge 

In the GLC gauge the averaging integrals introduced in Sect. 2 greatly simplify, especially 
in the case where the reference hypersurface T,(A) defines a geodesic observer. We will 
concentrate on the integrals (2.9) and (2.7) which, setting V = w,Vq = wq, are now given 
by: 



I(S; w , Ao; -) = J d 2 ddwdT^r\ \d T A\ 5(w - w )5(A - A )S , (4.8) 

7(5; w ;Ao) = f d 2 9dwd T J\^\ -j==5{w - wo)Q{A - A )S. (4.9) 

For a geodesic reference observer, with A = r (and Ao = To), we obtain: 

I (S; wo, to;-) = J d 2 9^(w ,T ,9 a )\ S{w ,T ,e a ), (4.10) 

I(S;wo;t ) = J d 2 edT^h(w ,T,e*)\ 6(r - T )S(w ,T,d a ). (4.11) 

The generalized Buchert-Elhers commutation rules introduced in Sect. 3 also simplify as 
illustrated below. 



Average on the 2-sphere embedded in the light-cone 
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Let us consider separately, as before, derivatives with respect to Aq and to Vq = up- 
starting with Aq, Eq. (3.4) can now be written as 

W*^ ) W„ oA ,(4.12) 

and for the special case where T,(A) defines a geodesic observer (with A = r and ^4o = 7"o)> 
it reduces to 

dr (S) WihT0 = {d T S) W0:T0 + (Sd r In - ld T In (S)^ . (4.13) 

\ / W ,T \ I W ,T 

Note that, in this last case with A = r, the partial derivative with respect to Aq reduces to 
a standard derivative with respect to the proper (cosmic) time if we consider the limit of a 
homogeneous FLRW Universe (see Eq. (4.4)). 

In the geodesic case we have also the following simplification for the derivative with 
respect to wq as given in Eq. (3.6): 

d W0 (S) W0 , T0 = (d w S + U a d a S) WOtTo + {Sd w ln ] /\ 1 \ 

w , T 



+ s 



d a u a + U a d a \nd\ 1 \ 



-(^ln^) (S) 



W ,T 

WQ, TO \ / Wo, TQ 



-( da U a + U a d a ln x /\ 7 \) (S) W(hTQ . (4.14) 

wo, TO 



Average on the truncated light-cone 

For a generic hypersurface we cannot specify {dA) 2 , and Eq. (3.5) just simplifies as 

* /qS A /d T S\ A ° I d T in^\-d w (dA-dA) \ Ao 

9A ° {S)w ° - \drA/ W o + \ drA / 



wo 



d T ln^\-d w (dA-dA)\ Ao 



d T A / 



(S)2 ■ (4.15) 



Wo 

WQ 



However, in the case where defines a geodesic observer (with A = r and Aq = tq) one 

obtains: 

dr (S)Z = (d T S) T ° o + (Sd r In ^\) T ° - (d T In (S)^ . (4.16) 

\ / WQ \ ' Wo 

Comparing with Eq. (4.13) we find, for this last case, the same commutation rule indepen- 
dently of the averaging prescription used. 

Finally, for a hypersurface E(.A) associated to a geodesic observer, we also have the 
following simplification of Eq. (3.7): 

dw (S)Z = (d w S + U a d a S)^ + (sd w ln^\) T ° 

\ / Wo 
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+ s 



d a U a + U a d a lnJ\ 7 \ 



WQ 



-{d w lnJ\ 7 \) (S) 



w 



\T0 
I WQ 



-/ da U a + U a d a lnJW\) T ° <S>£ . 

WQ 



(4.17) 



Hence, we obtain the same commutation rule (see Eq. (4.14)) even for derivatives with 
respect to wq. 



5 Some physical applications 

Information about the large scale structure of our Universe reaches us travelling along the 
null geodesies of a possibly inhomogeneous spacetime. Hence, an averaged description of the 
cosmological geometry should unavoidably make use of some light-cone averaging procedure, 
like those illustrated in the previous sections. In this section we will suggest two possible 
applications of the previous formalism: we will consider, in particular, the averaging of the 
luminosity-distance redshift relation and the averaging of the redshift drift parameter for a 
generic inhomogeneous cosmology. Note that in this section we will be dealing with the r 
coordinate both at the source and at the observer's position. They play the same role as 
Aq = to of the previous sections, and will be denoted by t s and tq respectively. 



5.1 The redshift to luminosity-distance relation 

It is well known that the scalar redshift parameter z, for a photon with momentum k^ 
emitted by a source S and received by the observer O, depends on the scalar product 
between k^ and the four-velocity n M of a local reference observer, and is given in general 
(with self-explanatory notation) by: 

Note that this object may be regarded as a bi-scalar, as it depends on the ratio of the same 
scalar quantity calculated in two different spacetime points. In the GLC gauge of Sect. 4, 
and in the physically interesting case of a geodesic observer (that we will consider hereafter) , 
we can use Eq. (4.6) and we thus obtain: 

1 + * = ^, (5-2) 

1 s 

where T = T(w , r , 0g) and T s = T(w ,t s ,6%). 

In a cosmological context we can define at least three different distances of physical 
interest (see e.g. [20]): d s , namely the angular distance of the source as seen from the 
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observer (in the FLRW limit of Eq. (4.4) we have d s = ra(t s ), where r is the comoving radial 
distance of the source); d , namely the angular distance of the observer "as seen from" the 
source (in the FLRW limit we have d Q = ra(to)); and di, the so-called luminosity distance. 
These three observational distances are always related to each other, independently of the 
given cosmological model, by redshift factors [15, 21]: 



as a consequence of the Etherington reciprocity law [22] . 

On the other hand, in a general spacetime, we can also identify (point-by-point on the 
2-sphere) the square of the angular distance of the source with the ratio of the intrinsic 
cross-sectional area element to the subtended solid angle element, according to d 2 = dA/dQ. 
(see e.g. [14]). In our GLC gauge we have = (0, — 1/Y,0,0), and, as already mentioned, 
the 2-sphere embedded in the light-cone (corresponding to the integration region of Eq. 
(2.9)) is orthogonal to the photon momentum, so that the cross-sectional area element is 
proportional to yjYy\d9 1 d9 2 . Therefore, d 2 s can be written as 



where p is the proper radius of an infinitesimal sphere centered around the observer, w = wq 
defines the past light-cone connecting source and observer, and r = t s defines the spacelike 
hypersurface normal to n^ 1 at the source position. Eq. (5.4) easily reduces to d s = ra(t s ) in 
the FLRW case of Eq. (4.4). 

It is appropriate to mention that in the present discussion we are neglecting the possible 
occurrence of caustics in the past light cone of our observer, which would affect the assumed 
relation between area distance and angular size distance. The possible relevance of this 
interesting effect, and its dependence on gravitational lensing and on the distribution of 
inhomogeneities has been addressed in detail in [23]. 

The relevant variables related to supernovae observations are the luminosity distance di 
and the redshift z. Their relation has been studied within a gauge invariant approach, for a 
linearly perturbed FLRW metric, in [24]. Averaging their values on the two-sphere embed- 
ded in the light-cone, and using the above results, we obtain, for a general inhomogeneous 
metric background, 



d L = (l + z)d = (1 + zfd. 



'Si 



(5.3) 




(5.4) 




jd 2 e^/\ 7 (w ,T S ,9-)\ [r 2 (w ,T ,e a )/T 2 (w ,T S ,9 a )] d s (w ,r s ,e a ) 

Jd 2 9^h(w ,r s ,9-)\ 



(5.5) 



where d s is given by Eq.(5.4), and 




J d 2 9^(w , t s , 9-)\ [T^p, r , 9 a )/T(w , t s , 9 a )} 
Jd 2 9^\ 7 (w ,r s ,9-)\ 



(5.6) 
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Since our coordinates should not be pathological near w = Wq and r = To we should have 
lim r _>. T0 T(w = wo,T,9 a ) = T(w = wq,tq) = To indepent of 9 a . Hence, in Eq. (5.5) and 
(5.6), the factor T(wo,TQ,O a ) = To behaves like a constant, leaving a well defined integral 
of a scalar object over the 2-sphere. 

The above expressions can be used to see how the usual redshift to luminosity-distance 
relation gets affected by the inhomogeneities of the cosmological geometry. To this purpose 
it is also possible (and probably more convenient) to consider averages of the luminosity 
distance on constant redshift surfaces 3 . In that case the geodesic observers measuring the 
redshift do not coincide any longer with the "observers" associated to the flow lines of 
the reference hypersurface Ti(A), chosen to specify the averaging region. Even in such a 
case, however, our general formalism provides a rather simple averaging prescription by 
identifying the scalar A in Eq. (2.12) with k^u^ (where now is the velocity of a geodesic 
observer), and the quantity to be averaged with oIl (or with some convenient observable 
related to it). 

The straightforward calculation simplifies considerably in our gauge giving, for instance, 

, , x n , , 2 Jd 2 e^\ 1 (wo,r(z,w ,e a ),e a )\d s {wo,T(z,w ,e a ),e a ) 

{dL)w ,z = {1 + Z) n— aa \ =dET\ ' 

where T(z,wo,9 a ) is the solution of: 

T(ti»,T,9») 1 



To 1 + ,- < 5 ' 8 > 

We plan to come back to an explicit application of our formalism to this problem in the 
near future. 

5.2 Redshift drift 

The redshift drift (RSD) is defined as the rate of change of the redshift of a given source with 
respect to the observer's proper time. Since both the observer and the source simultaneously 
evolve in time, the relevant hypersurfaces for the RSD effect will consist of two light-cones 
with two different bases (see Fig. 2). Assuming that source and observer are in geodesic 
motion and have negligible peculiar velocity, the RSD in a FLRW Universe, and in the 
proper-time interval Ato, can be simply expressed as (see e.g. [13]): 

^ = (l + z)H -H s = a -^, (5.9) 
Ato a s 

where Hq is the present value of the Hubble parameter while H s is its value at the emission 

time. Clearly the RSD effect is a good indicator of cosmic acceleration as a function of z 

which does not use hypothetical standard candles. 



3 We would like to thank Misao Sasaki for stressing this point with us. 
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Figure 2: A graphic illustration of the redshift drift effect. A possible variation of the cos- 
mological expansion rate is detected by comparing observations performed on two different 
past light-cones. 



In order to generalize this expression to a general inhomogeneous Universe (and attempt 
some averaging of it) let us consider again the geodesic observer of the GLC gauge, with 
coordinates = (w,T,6 a ). As shown in Eq. (5.2), z is in principle a function of seven 
independent variables, namely of wo, tq, 9q, t s and 6® (note that, since we are assuming 
that source and observer are on the same light-cone w = wq at r = tq, they will be both 
on the light-cone w = wq + Awo at t = ro + A to, see Fig. 2). We can also note that To 
is independent of 9q (see the comment after Eq.(5.6)). As a consequence we have only five 
independent contributions to the variation of 1 + z, and we can write, in general: 



. d(l + z) A 8(1 + z). 0(1 + *) a 5 ( 1 + *)a/. 



(5.10) 



As shown in Eq.(4.4), in the homogeneous limit the coordinate r goes to the proper (cos- 
mic) time t of the synchronous gauge. So, locally around our geodesic observer, we choose 
to evaluate the redshift drift Az with respect to his/her time parameter r, and we need, 
to this purpose, an explicit relation between Atq and the variation of the other coordinates 
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involved. For a geodesic observer with = —d^r we have x M ~ rc/ 4 = (1/T, 1, U a /T), so 
that we can express Ax' 1 in terms of Aw, at all times, as: 



At = TAw 



A6 a = U a Aw. 



(5.11) 



Using Eq. (5.11) we find Awq = Ato/To = Ar s /Y s , and the final result for the RSD can 
be written in terms of z and its derivatives as: 



where: 



1 dTn 



(5.12) 



(5.13) 



T 3t 

(let us note again that d WQ acts on both metric factors T contained in (1 + z)). This result 
is valid for a general inhomogeneous metric, and can be compared, as a useful consistency 
check, with a similar result previously obtained in the particular case of a spherically sym- 
metric geometry [12]. If we move from our coordinates to the adapted coordinates used 
in [12] we find that our expression (5.12) exactly reduces to the expression for the RSD 
reported in Eq. (5) of Ref. [12]. 

The quantities appearing in Eq. (5.12) can now be averaged over the 2-sphere embedded, 
at r = t s , in the light-cone w = wq. Using our prescription based on Eq. (4.10) we find 
that both To and Hq can be taken out from the averaging integrals (which are performed 
over the variables 9 a ), and we obtain: 

(Az) W0 , Ts ={1 + z) h q + _L {dw{1 + z))w(hTs + {drHl + z))mhTs + _L (u a d a (l + z)) m , Ts . 

A-T~0 I i-0 

(5.14) 

For the derivatives performed with respect to the hypersurface parameters w and r we can 
now apply the Buchert-Ehlers commutation rules in the simplified form of Sect. 4.2. Using 
in particular Eqs. (4.13) and (4.14) we are lead to: 



(A*) 



At 



(1 + z)w ,t s H + ^d W0 (l + z) WOjTs + d Ts (ln(l + z)) WOjTs 
lo 

+Qw(z) + Qt{z), 



where 



Q T {z) = - { ln(l + z)d T In ^/| 7 | ) + (ln(l + z))^ ( d T In ^/| 7 | 

WQ,T S \ / W0,T s 



(5.15) 



(5.16) 



Qw(z) 



1 



- <(l + z)d w \nJ\ 7 \ 



T x 

+ (1 + z) WOtTs 



WQ,T S 



(l + z) 



d a U a + U a d a ln^\ 



WQ, T S 



d w hiy0\ +(d a U a + U a d a ln^, 



wo,r s 



(5.17) 
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These last terms, together with the other terms of Eq. (5.15) control the difference be- 
tween the averaged redshift drift and the drift for the averaged redshift and represent a 
light-cone analog of the backreaction terms due to the geometric inhomogeneities, first 
computed by Buchert [5] in the context of three-dimensional spacelike averages. Unlike in 
that case, however, our backreaction is physically controlled by the geometric dynamics of 
a two-dimensional surface, with metric ^ a b- Also, the possible physical meaning of these 
backreaction terms is in principle strongly model-dependent, and their physical effects are 
to be explicitly computed for any given model chosen to parametrize deviations from FLRW 
geometry. 

For a homogeneous FLRW metric the averages disappear, and the backreaction terms 
Q T and Q w are identically vanishing. In that case Hq = H$ and, using the homogeneous 
limit of Eq.(4.4), one also finds d w (l + z) = and <9 T ln(l + z) s = —H s , so that Eq. 
(5.15) gives back the result (5.9). For a general inhomogeneous model, however, several 
new contributions appear. It would be very interesting to see whether such backreaction 
effects are relevant for determining the kinematics of the Universe (and the equation of state 
of its energy components as a function of z) from the forthcoming RSD experiments (see 
e.g. [25]). We plan to come back to this question elsewhere [26], at least within specific 
inhomogeneous models such as the Lemaitre-Tolman-Bondi (LTB) Universe [27]. 

6 Conclusions 

In this paper we have made a first step towards defining suitable covariant and gauge 
invariant light-cone averages that should be relevant to analyze the effects of inhomogeneities 
on astrophysical observables related to light-like (massless) signals. We were able to define 
such averages both on light-like hypersurfaces and on ordinary 2-surfaces embedded on a 
specific light-cone (e.g. our past one). We have also written down the generalized version 
of the Buchert-Ehlers commutation rules (between averaging and differentiation) that hold 
for all such light-cone averages. 

One obvious problem to which we would like to apply our technique is that of cosmic ac- 
celeration. This is actually nothing but a measured relation between redshift and luminosity 
distance which can only be explained in a FLRW cosmology if the Universe underwent a 
late-time accelerated expansion. Within General Relativity it implies, in turn, the existence 
of a large dark energy component in the cosmic fluid. Our formalism allows in principle to 
study directly the luminosity- redshift relation for a given non homogeneous Universe (e.g 
of the LTB type [27]) or, even more interestingly, for a Universe containing a stochastic 
spectrum of inhomogeneities like those that originate from inflationary cosmology. 
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In order to prepare the ground for such investigations we have introduced some adapted 
"geodesic null" coordinates that allow to express in a simple way redshift and luminosity 
distance in terms of local metric components (i.e. calculated at the source's position). Fur- 
thermore, our expressions for both the averages and for their derivatives take a considerably 
simpler form in these coordinates. 

We have also applied our formalism to the case of the so-called redshift drift, a quantity 
which, without any need of using standard candles, should be sensitive to a possible ac- 
celeration of the scale factor at different redshift values, and thus could distinguish among 
different forms of dark energy. Once more, however, one should address the question of 
how inhomogeneities could modify the FLRW relation between redshift drift and expansion 
rate. By using the geodesic light-cone gauge we were able to find the formal modification 
of the relation and show that it contains, among others, Buchert-Ehlers-like commutator 
terms. The physical interpretation of such terms is, however, model-dependent, and their 
possible effect are to be extracted for any given model with an explicit computation. 

Unlike in the case of spacelike hypersurfaces we have found that the most useful averages 
for physical applications are related to surface averages. However, it is easy to see that the 
quantities to be averaged are themselves sensitive to inhomogeneities over the whole past 
light-cone hypersurface even when they can be expressed in terms of quantities living on 
the surface. So far we have not found physical examples in which a true average over the 
full light-cone is involved, and we strongly suspect that there isn't one. 

In order to find out whether the formalism we have introduced is actually useful one will 
have to perform detailed calculations within realistic and explicit models for the relevant 
inhomogeneities. These can be either abstract mathematical ones, like LTB, swiss-cheese, or 
fractal models, or realistic ones based on our present knowledge of the large scale structure. 

Another interesting exercise would be to compute the magnitude of several effects to 
second order in a perturbed FLRW Universe. Even if, as most people suspect, the effect 
of inhomogeneities will turn out to be small and not to change the conclusions that follow 
from simple FLRW consideration, the result of such an investigation will be important in 
that it will sharpen considerably the conclusion that dark energy is indeed unavoidable. 
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